Antiproton and Positron Signal Enhancement in Dark Matter Mini-Spikes Scenarios 
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The annihilation of dark matter (DM) in the Galaxy could produce specific imprints on the spectra 
of antimatter species in Galactic cosmic rays, which could be detected by upcoming experiments 
such as PAMELA and AMS02. Recent studies show that the presence of substructures can enhance 
the annihilation signal by a "boost factor" that not only depends on energy, but that is intrinsically 
a statistical property of the distribution of DM substructures inside the Milky Way. We investigate 
a scenario in which substructures consist of ~ 100 "mini-spikes" around intermediate- mass black 
holes. Focusing on primary positrons and antiprotons, we find large boost factors, up to a few 
thousand, that exhibit a large variance at high energy in the case of positrons and at low energy 
in the case of antiprotons. As a consequence, an estimate of the DM particle mass based on the 
observed cut-off in the positron spectrum could lead to a substantial underestimate of its actual 
value. 
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Many current and upcoming experiments are aimed 
at detecting the products of annihilation of dark matter 
particles in the Milky Way halo, or in external galaxies 
- see Refs. [1, 2, 3] for recent reviews. In fact, the iden- 
tification of an annihilation signal would not only pro- 
vide the proof of the particle nature of DM, but it would 
also open a new window to Physics beyond the Standard 
Model. Among indirect searches, particular attention has 
been devoted to the case of photons, since they propa- 
gate along straight lines, and they are not significantly 
absorbed within the galaxy at energies of order 100 GeV 
and below. The prospects for detecting gamma-rays from 
DM annihilation with the space telescope GLAST [4] and 
with Air Cherenkov telescopes such as CANGAROO [ ■".], 
HESS [6], MAGIC [7] and VERITAS [8], have been ex- 
tensively discussed in literature - e.g. [52] and references 
therein. Alternatively, one may search for neutrino fluxes 
from DM captured at the center of the Sun or of the 
Earth - see reviews above for further details - with large 
neutrino telescopes such as IceCube [9] and Antares [10]. 

One can also focus on other final states of DM an- 
nihilation, such as antimatter particles, in particular 
positrons [11, 12, 13, 14, 15, 16, 17, 18, 19] and antipro- 
tons [20, 21, 22, 23]. Unlike photons and neutrinos, these 
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are charged particles that are sensitive to the galactic 
magnetic field and suffer energy losses. Their propaga- 
tion must be described by diffusion, so that the magnetic 
properties of the Milky Way must be known, which is not 
completely true yet. Nevertheless, it is possible to make 
interesting predictions, which are now particularly rele- 
vant, as the recently launched satellite PAMELA [24] and 
the forthcoming AMS02 [25] experiment are expected to 
provide new data and test a significant portion of the 
parameter space of DM particles. 

In most cases, the prospects for the detection of anti- 
matter fluxes are not particularly sensitive to the overall 
shape of the DM profile, but they do depend strongly on 
the "dumpiness" of the profile. In fact, the annihilation 
rate being proportional to the square of the number den- 
sity, the presence of small-scale inhomogeneities leads to 
larger fluxes for the same value of the average DM den- 
sity. This is often parameterized by a "boost factor" , 
that will be defined more precisely below. 

Recenty, a new scenario has been discussed in the 
literature, where the formation of Intermediate Mass 
Black Holes, i.e. Black Holes with mass M in the range 
10 2 < M/Mq < 10 6 , leads to the formation of DM over- 
densities, called "mini-spikes" , that might be observed as 
point sources of gamma-rays [40] and neutrinos [26]. Al- 
though somewhat speculative, the scenario has the undis- 
puted virtue of making specific predictions on the number 
and luminosity of these objects, that could be observed 
or ruled out in the near future with the upcoming gener- 
ation of space and ground based experiments. Here, we 
compute the boost factor in mini-spikes scenarios, and 
argue that they would lead to a dramatic enhancement 
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of antimatter fluxes, bringing them within the reach of 
current and upcoming experiments. 

The paper is organized as follows : in Section I, we re- 
view the motivations for the mini-spike scenario and de- 
rive their density profile and normalization. Section II is 
devoted to a general discussion of the antimatter cosmic 
ray signal generated by mini-spikes. The average value 
and the variance of the boost factor at the Earth are de- 
fined and special attention is paid to the inner structure 
and galactic distribution of the mini-spikes. Positrons 
and antiprotons are respectively addressed in Sections III 
and IV. Results from Monte-Carlo simulations are pre- 
sented and discussed in the light of analytic calculations. 
The energy dependence of the boost variance turns out 
to be quite different for positrons and antiprotons. These 
species have very different modes of propagation inside 
the Milky Way. A few realistic models for the DM par- 
ticles are implemented in Section V. The positron and 
antiproton signals which they yield are featured and il- 
lustrate the typical antimatter signatures to be expected 
should mini-spikes populate the galactic halo. We finally 
conclude in Section VI and make a few suggestions for 
future developments. 

I. MINI-SPIKES AS DM SUBSTRUCTURES 

Black Holes (BHs) can be broadly divided in three dif- 
ferent classes - see e.g. [27] for a review : 

• Stellar Mass BHs, with mass M < 100 M Q . They 
are typically remnants of the collapse of massive 
stars. Recent simulations suggest that the upper 
limit on the mass of these objects is as low as 
~ 20 M Q [28]. There is robust evidence for the 
existence of these objects, coming from the obser- 
vation of binary systems with compact members 
whose mass exceeds the critical mass of Neutron 
Stars. For a review of the topic and the discussion 
of the status of the observational evidence for Stel- 
lar Mass BHs sec e.g. [29] and references therein. 

• Supermassive BHs (SMBHs). with mass M > 
10 6 M Q lying at the centers of galaxies, including 
our own. Their existence is also well-established - 
see e.g. Ref. [30] - and intriguing correlations are 
observed between the SMBHs mass and the proper- 
ties of their host galaxies and halos [31, 32, 33, 34, 
35, 45]. From a theoretical point of view, a pop- 
ulation of massive seed black holes could help to 
explain the origin of SMBHs. In fact, observations 
of quasars at redshift z ss 6 in the Sloan Digital sur- 
vey [36, 37, 38] suggest that SMBHs were already 
in place when the Universe was only ~1 Gyr old, 
a circumstance that can be understood in terms of 
a rapid growth starting from "massive" seeds - see 
e.g. Ref. [39]. This leads to the third category: 

• Intermediate Mass BHs (IMBHs), with mass 10 2 < 
M/Mq < 10 6 . Scenarios that seek to explain 



the properties of the observed supermassive black 
holes population result in the prediction of a large 
population of wandering Intermediate Mass BHs 
(IMBHs). Here, following Ref. [40], we consider 
two different formation scenarios for IMBHs. In 
the first scheme (A), IMBHs form in rare, over- 
dense regions at high redshift, z ~ 20, as rem- 
nants of Population III stars, and have a charac- 
teristic mass-scale of a few 10 2 M [41] - a simi- 
lar model was investigated in Ref. [42, 43, 44]. In 
this scenario, these black holes serve as the seeds 
for the growth of supermassive black holes found 
in galactic spheriods [30]. In the second scenario 
(B), IMBHs form directly out of cold gas in early- 
forming halos and are typified by a larger mass 
scale of order 10 5 M Q [45]. During the virialization 
and collapse of the first halos, gas cools, collapses, 
and forms pressure-supported disks at the centers 
of halos that are massive enough to contain a large 
amount of molecular hydrogen. In halos which do 
not experience any major mergers over a dynami- 
cal time, a protogalactic disk forms and can evolve 
uninterrupted. At this stage, an effective viscosity 
due to local gravitational instabilities in the disk 
leads to an inward mass transfer and outward an- 
gular momentum transfer, until supcrnovae in the 
first generation of stars heat the disk and terminate 
this process. By the time the process terminates, a 
baryonic mass of order ~ 10 5 Af Q loses its angular 
momentum and is transferred to the center of the 
halo, leading to the formation of an object that may 
be briefly pressure-supported, but that eventually 
collapses to form a black hole. In Fig. 4 wc show 
the distribution of IMBHs in the latter scenario, as 
obtained in Ref. [26]. 

The effect of the formation of a central object on the 
surrounding distribution of matter has been investigated 
in Refs. [46, 47, 48, 49] and for the first time in the 
framework of DM annihilations in Ref. [ ] . It was shown 
that the adiabatic growth of a massive object at the center 
of a power-law distribution of DM, with index 7, induces 
a redistribution of matter into a new power-law (dubbed 
"spike") with index 

7 sp = (9-2 7 )/(4- 7 ) . (1) 

This formula is valid over a region of size R sp w 0.2 rBH , 
where tbh is the radius of gravitational influence of the 
black hole, defined implicitly as M(< tbh) = Mbh, 
where M(< r) denotes the mass of the DM distribution 
within a sphere of radius r, and where Mbh is the mass 
of the Black Hole [51]. The process of adiabatic growth 
is in particular valid for the SMBH at the galactic cen- 
ter. A critical assessment of the formation and survival 
of the central spike, over cosmological timescales, is pre- 
sented in Refs. [52, 53] - see also references therein. Adi- 
abatic spikes are rather fragile structures, that require 
fine-tuned conditions to form at the center of galactic 
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halos [54], and that can be easily destroyed by dynami- 
cal processes such as major mergers [55] and gravitational 
scattering off stars [53, 56]. 

It was recently shown that a p ot r~ 3 / 2 DM overden- 
sity can be predicted in any halo at the center of any 
galaxy old enough to have grown a power-law density 
cusp in the stars via the Bahcall-Wolf mechanism [57]. 
Collisional generation of these DM "crests" - Collision- 
ally REgenerated STructures - was demonstrated even 
in the extreme case where the DM density was lowered 
by slingshot ejection from a binary supermassive black 
hole. However, the enhancement of the annihilation sig- 
nal from a DM crest is typically much smaller than for 
adiabatic spikes [57]. 

Here we focus our attention on mini-spikes around 
IMBHs, and we recall their properties, following closely 
Ref. [40]. The "initial" DM mini- halo - that is, the DM 
distribution prior to black hole formation - can be well 
approximated with a Navarro, Frenk, and White (NFW) 
profile [58] 

p(r)=Po(—) U + 2 ■ (2) 

\r S cJ \ r sc J 

The normalization constant po- and the scale radius r sc , 
can be expressed in terms of the virial mass M v { r of the 
halo at the time when the IMBH formed, and of the virial 
concentration parameter c V i r _ see Ref. [40] for futhcr 
details. Alternatively, we could have chosen the more 
recent parameterization proposed by Navarro et al. [59] 
- see also Refs. [60, 61]. However, this profile implies 
modifications at scales smaller than those we are inter- 
ested in, where the profile is anyway modified by the pres- 
ence of the IMBH. We assume that the black holes form 
over a timescale long enough to guarantee adiabaticity, 
but short compared to the cosmological evolution of the 
host halo. In fact, the condition of "adiabaticity", funda- 
mental to grow mini-spikes, requires that the formation 
time of the black hole is much larger than the dynamical 
timescale at a distance tbh from the black hole, where 
Tbh — GMbh/o 2 is the radius of the sphere of gravita- 
tional influence of the black hole, and a is the velocity 
dispersion of DM particles at tbh- In practice, Tbh is 
estimated by solving the implicit equation 

rrBH 

M(<r BH )= / Ptot (r)4irr 2 dr = 2M BH ■ (3) 
Jo 

In the following, we only consider the formation scenario 
B, which leads to higher values for the boost factor. For 
a representative case in scenario B, with Mbh = 10 5 M Q 
and M v ir,/ = 1O 8 M , this gives r B H/r BC ~ 0.04. In 
scenario B, the black hole formation time is set by the 
timescale for viscous angular momentum loss and is lim- 
ited by the evolutionary timescale of the first stars and 
the gravitational infall time across the gaseous disk, 
which is of order Myr - see Ref. [62] for a detailed discus- 
sion of timescales. The relevant timescale for the mass 
build up of the IMBH is then t ev ~ 1 - 20 Myr. 



II. THE ANTIMATTER COSMIC RAY SIGNAL 
PRODUCED BY MINI-SPIKES 

A. Fluxes at the Earth 

The self-annihilations of the DM particles \ concealed 
in the Milky Way halo produce positrons and antipro- 
tons, a small fraction of which may reach the Earth. 
The local production rate of these antimatter cosmic ray 
species will be generically denoted by 

P(x) = 5 (a ann v) J^Mj g (E s ) AE S . (4) 

The factor S is equal to 1/2 if the DM particles \ are 
Majorana fermions whereas it is 1/4 in the case of Dirac 
species, provided that a perfect matter-antimatter sym- 
metry between x an d X holds in that case. Until Sec- 
tion V, we will focus our discussion on monoenergetic 
positrons and antiprotons with energy E$ defined up to 
an energy bandwidth of IS.Es- In Kaluza-Klein theories 
of dark matter particles, the branching ratio into an e + e~ 
pair can be as large as B e ± ~ 20 % [63] and the energy 
distribution g{Es) of Eq. (4) is equal to B e ± S(Es — m x ), 
where m x stands for the mass of the DM species. In the 
case of antiprotons, energy losses may be neglected and 
we need only to consider a specific value of the initial 
energy E$- A continuous injection spectrum will be con- 
sidered in Section V, as is needed for realistic models of 
antimatter production. 

The probability for a cosmic ray species injected at 
x to propagate and reach the Earth with energy E is 
described by the Green function 

G(x)=G{E,Q^E s ,x) . (5) 

The propagator G(x) has been thouroughly discussed 
in the literature - see Refs. [16, 19] for positrons and 
Refs. [64, 65] for antiprotons. The annihilations of DM 
particles yield the following flux at the Earth 

^{E)=S( G(x) (^M) 2 d 3 x , (6) 

J DM halo I PO ) 

which is expressed in units of cm" 2 s~ x sr _1 GcV -1 . 
The factor S encodes informations on the specific model 
selected to describe the DM species x such, as its mass m x 
and total annihilation cross section (<T ann v). It is defined 

as 

S = /■ /3(E) (a ann v) { - ) g(E s )AE s . (7) 
4tt lm x ) 

The velocity of the cosmic ray particles is denoted by 
P(E) and po is a reference density that can be chosen 
at will, eventually disappearing from the final result. It 
can be given a physical meaning, as we will see when we 
define the annihilation volume £. 
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B. The framework of the statistical analysis 

In the absence of any substructure, the DM distribu- 
tion p s is smooth and yields a cosmic ray flux 
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(8) 



In the IMBHs scenario, the DM distribution is given by 
the superposition p = p s + Sp where a new component 
Sp ^> p s accounts now for the DM trapped inside the 
mini-spikes. That component contributes for a fraction 
of ~ 10~ 5 to the Galactic dark matter. As the typical 
distance over which positrons and antiprotons propagate 
in the Milky Way is much larger than the average size of 
the mini-spikes, the cosmic ray flux - generated by the 
total density p - simplifies into the sum 



<t>(E) 




(9) 



where the contribution from the i-th object is ifi = 
S x G(xj) x £j. That substructure produces as many 
positrons and antiprotons as if the entire volume 
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i-th mini-spike 



Po 



d 6 x 



(10) 



were filled with the constant DM density po- The mini- 
spike mass Mi and intrinsic boost Bi are related to the 
annihilation volume £j through 
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Bj Mj 
Po 



(11) 



Note that the exotic flux tpi has no well-defined depen- 
dance on the mass Mj and the intrinsic boost Bi if these 
quantities are considered individually, since (fi only de- 
pends on the annihilation volume £j. 

Should we know the exact location and annihilation 
volume of each mini-spike, we would unambiguously de- 
rive the cosmic ray flux 4>(E). This is not the case, 
and for that matter wc do not even know the number 
Abh of these objects. Very different distributions of the 
mini-spike population throughout the galactic halo are 
possible. The positron and antiproton spectra at the 
Earth are therefore affected by an uncertainty (which 
could be called a cosmic variance) related to this lack 
of knowledge. A whole set of halo realizations has to be 
considered, and specific semi-analytical tools have been 
developped [19] to address this issue. The boost fac- 
tor B = 4>/4> s at the Earth is not unique and must be 
treated as a random variable. To do so, we have built a 
Monte-Carlo simulation nurtured by ~ 200 different re- 
alizations of the mini-spike population which have been 
obtained in Ref. [40] by evolving an initial distribution 
of IMBHs orbiting in the Milky Way halo and by allow- 
ing the associated DM mini-halos to be tidally destroyed 



during close encounters. The distribution of the num- 
ber Abh of surviving mini-spikes is presented in Fig. 1. 
In our Monte-Carlo simulation, Abh is randomly drawn 
according to a gaussian distribution. About 100 objects 
are expected within a 100 kpc galactocentric radius and 
~60 of them populate the Galactic diffusive halo. 
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RMS 20.73 
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FIG. 1: Distribution of the Monte-Carlo realizations of the 
galactic mini-spike population - extracted from Ref. [40] - as a 
function of the number JVbh of objects within a galactocentric 
radius of 100 kpc. 

If we set Abh equal to its mean value, we may derive an- 
alytically the average value and the variance of the boost 
factor B according to the method presented in Ref. [19]. 
We did not find any significant correlation between the 
position x and annihilation volume £ of the mini-spikes 
in the simulations of Ref. [40], so that they seem to be 
distributed independently of the others. We readily infer 
an average boost factor of 
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The integral X stands for the convolution of the propa- 
gator G with the smooth DM distribution p s 
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In the following, a NFW profile is assumed, where a 
value po — 0.3 GeV/cm 3 is considered for the Solar dark 
matter density. The mini-spike annihilation volume and 
galactic position are respectively distributed according 
to the probability functions q(£) and p(x) which are pre- 
sented in Sections II C and II D. These probability laws 
allow to define the average values 



(G n ) = 



{G(x)}™ p(x) d 3 x 



(14) 
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(15) 



In practice, we will only be concerned with n = 1 and 2. 
The boost scatter <tb may be derived from the relation 
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(16) 



where the variance oy of the random flux component <j) r 
is given by 
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(17) 

where cr^, ctg and cr/v stand for the variances of the re- 
spective distributions. In Sections III and IV, we will 
check the consistency of our Monte-Carlo results in the 
light of the analytic expressions (12) and (16). 



C. The inner structure of mini-spikes 

In the case where the DM profile before the formation 
of the IMBH follows the commonly adopted NFW distri- 
bution [58], the final DM density p(r) around the IMBH 
will be described by a power law r~ 7 / 3 in a region of size 

Rsp' 




\og{R) 



FIG. 2: Schematic representation of the inner structure of a 
mini-spike 

At larger distances, the DM distribution has not been 
modified by the accretion onto the IMBH and the density 
still falls down as r . This envelope does not contribute 
significantly to the DM annihilation and its associated 
production of positrons and antiprotons. On the con- 
trary, the DM density steeply increases below R sp and 
annihilations themselves set an upper limit to it 
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where r is the time elapsed since the formation of the 
mini-spike. We denote by R c the "cut-off" radius below 
which the mini-spike core extends with uniform density 
Pmax- For a typical value of r = 10 Gyr, we find 



Pmax = 8.752 x 10 1 



(19) 
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The relative extension of the mini-spike mantle (R c < 
r < Rsp) with respect to the core (r < R c ) is given by 
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where p sp = p(R sp ) is the DM density at the surface 
r = R sp of the mantle. Integrating Eq. (10) over the 
inner structure of the mini-spike leads to the annihilation 
volume 
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(21) 



Due to the different mini-spike features (much less steep 
DM profile), the case of scenario A leads to much lower 
annihilation volumes. This case is no more mentioned in 
the following. 
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FIG. 3: The probability law q(£) for the annihilation volume 
has been derived from the Monte-Carlo results of Ref. [40]. 
We found no correlation with the mini-spike position. 

We have derived the probability function for £ from the 
Monte-Carlo results of Ref. [40] where the mini-spike ra- 
dius R sp and external density p sp are provided for each 
object. This distribution - featured in Fig. 3 - is charac- 
terized by very large values of the annihilation volume, 
with a tail extending up to £ ~ 1.2 x 10 7 kpc 3 . On aver- 
age, with a mini-spike radius R sp = 2.84 pc and density 
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FIG. 4: Radial distribution of the mini-spikes, as extracted 
from the numerical results of Ref. [40]. 

p sp = 48.51 M Q pc~ 3 , we infer a typical annihilation vol- 
ume £ ~ 3.3 x 10 6 kpc 3 for a benchmark cross section 
(camiu) — 3 x 10~ 26 cm 3 s _1 , a mass m x = 1 TcV and 
Po = Pq = 0.3 GcV cm" 3 . The distribution of £ hap- 
pens to be log-normal. In our Monte-Carlo program, we 
simulate a large number of different IMBH halo popula- 
tions. For each mini-spike, the annihilation volume £ is 
randomly chosen according to the distribution. 

Notice finally that the individual mini-spike contribu- 
tion ipi scales as the product iS x and is eventually 
proportional to 

^mini-spike oc (tT ann v) 2/7 rt\ 9/7 ■ ( 22 ) 

The cosmic ray signal depends weakly on the annihila- 
tion cross section since a decrease of the cross section is 
partially compensated by a higher annihilation volume. 
This issue is addressed with more details in Section V, 
where the influence of the injection spectrum g(Eg) will 
also be considered. 



D. The galactic distribution of mini-spikes 

The Monte-Carlo simulation of IMBHs formation and 
evolution is used to derive the distribution function p(r) 
for the galactocentric radius of each mini-spike that has 
survived tidal disruption. This function gives the proba- 
bility to find an object at galactocentric distance r within 
a spherical shell of thickness dr. It is expressed in units of 
kpc" 3 and normalized to unity within the inner 100 kpc. 
This volume distribution is presented in Fig. 4 together 
with the parameterized function used in the following to 
estimate the boost factors. 

The mini-spike number density drops quickly beyond 
a few hundred parsecs from the galactic center. This is 



because encounters and mergers are much more frequent 
in this region, and they are very efficient in disrupting 
dark matter spikes. For those which have survived, the 
logarithmic slope of the radial profile varies from ~ 1.8 
in the inner region to ~ 3.8 outward. This distribution 
is more peaked than for a typical NFW profile whose 
corresponding volume probability is also drawn in Fig. 4 
for comparison. 



III. RESULTS FOR POSITRONS 

This section is devoted to the positron signal. For 
pedagogical purposes, we first consider a fiducial case 
in which dark matter consists in m x = 1 TeV particles 
fully annihilating at rest in electron-positron lines. In 
this very simple frame, we scrutinize the positron flux 
enhancement that mini-spikes produce with respect to a 
single smooth distribution of dark matter in the Galaxy 

The propagator for positrons is given by [19]. The en- 
ergy at which a positron is detected on Earth depends on 
its propagation history. It is correlated to the distance 
over which it has diffused from the source. As a result, 
the flux at a given (detected) energy depends on the spa- 
tial distribution of the sources, i.e. on dumpiness. In 
this section, we study this effect in a more quantitative 
way. 



A. The Monte-Carlo simulation 

We have performed Monte-Carlo simulations of the 
mini-spike spatial distribution in the halo. Each halo con- 
figuration yielded a value for the positron flux 4> r = ^2 P* > 
from which the average boost factor and the correspond- 
ing variance could be computed. The probability density 
functions for the number Nbh of IMBHs, their position 
Xj and annihilation volume £j were based on cosmological 
simulations from Ref. [40]. The corresponding probabil- 
ity laws have been presented in Section II. The propa- 
gated flux is compared to the flux produced by a smooth 
NFW profile to determine the boost factor. 



B. The numerical results 

Monte-Carlo simulations involving 10 6 realizations of 
the IMBHs population have been performed. The mean 
value and the variance of the resulting boost factors (for 
positron flux) at the Earth is sketched in Fig. 5 as a 
function of energy. The shaded area show the region 
where the boost factor is expected to lie, with a 1-cr and 
2-er confidence level. 
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FIG. 5: Results from the Monte-Carlo simulations of the 
IMBHs population inside the Milky Way are compared to 
the analytical computations of the effective boost factor and 
its dispersion, for m x = 1 TeV. 

The yellow (grey) areas correspond to the 1 a region, 
the lighter one being obtained by fixing the annihilation 
volumes to their mean expected value, while the darker 
one corresponds to the general case for which mini-spikes 
have values of £j drawn according to their true probability 
distribution g»(£). In both cases, the dot-dashed curves 
stand for the 1 o contours obtained analytically. These 
curves are in good agreement with the ones obtained from 
the Monte Carlo, the small increase of the variance for the 
Monte Carlo with respect to the analytical expectation 
is due to the variation of the black holes number from 
one Monte Carlo realization to another, while this effect 
is not implemented in the analytical determination of the 
boost factor. The fact that these curves are fairly close 
to each other confirms that the dispersion of the number 
-/Vbh of IMBHs in the Milky Way poorly influences the 
final dispersion of the boost. This figure also shows that 
the boost factor can be very large, the expected value 
being of order 8000 (for a DM particle mass of 1 TeV). 

Fig. 5 shows that the variance increases with energy, 
and becomes very large for E > 20 GeV. It clearly ap- 
pears that the boost factor obeys two different statistical 
regimes. This is because the diffusive range of positrons 
depends on energy and a positron emitted with energy E s 
looses energy as it propagates outwards from the source. 
The typical propagation scale for a positron injected with 
energy E s and detected with energy E is given by 

Xd = (e 3 - 1 - el" 1 ) /(I - 8) . (23) 

where Kq and 5 are respectively the normalization and 
the logarithmic slope of the diffusion coefficient, and 
e = (E/l GeV). We have considered the median param- 
eters given in [21], K = 0.0112 kpc 2 Gyr" 1 , L = 4 kpc, 



V c = 12 km/s and S = 0.7. The typical timescale for 
energy loss of 1 GeV positrons is te ~ 10 16 s. For a 
1 TeV injected energy, we infer a propagation length 
\d — 6.9 kpc x -\/e -0 ' 3 — 0.12, which ranges from 0.1 
kpc at a detected energy of 990 GeV to 6.4 kpc at 1 
GeV. 

At high energy, this diffusive range is very short, and 
positrons will be detected only provided there happens to 
be a mini-spike very close to us. The probability of this 
event is low at energy close to E s , but the contribution 
to the flux is high if it happens. The trade-off between 
these two trends leads to a constant expectation value 
for the boost factor, but a large shotnoise-like variance. 
Conversely, at low energy the diffusive range is large, and 
the fluctuations in the number of sources contributing to 
the flux becomes relatively small. 

The distributions of the boost factor at 1 GeV ob- 
tained with Monte Carlo simulations are displayed in the 
left panel of Fig. 6. The red histogram corresponds to 
the case of identical mini-spikes with fixed annihilation 
volume. It can be compared to the analytical estimate of 
the boost factor distribution (dashed curve) in the limit- 
ing case where the number of IMBHs is large enough for 
the central-limit theorem to be valid. If it were so, the 
distribution of the boost factor would tend to be gaus- 
sian. One can see that it is almost the case for identical 
objects. Though unrealistic, the case in which ten times 
more IMBHs populate the Milky Way leads actually to 
a gaussian distribution for B(l GeV) (not shown here). 
The statistics obtained in the general case is plotted in 
black, one can see that the effect of having different an- 
nihilation volumes for different mini-spikes shifts the dis- 
tribution away from the gaussian behavior. 

The right panel of Fig. 6 displays the distribution of 
the boost factor at 300 GeV. Again, the red histogram 
corresponds to the identical IMBHs case and should be 
compared to the analytical estimate (dashed curve) . The 
latter includes only one IMBH inside the sensitivity vol- 
ume so that the extra events of the red curve correspond 
to the (very) rare situations in which two mini-spikes 
contribute to the signal. This feature is not present in 
the general case, for which the statistics is displayed in 
black : it is erased by the random choice of £. 



IV. RESULTS FOR ANTIPROTONS 

A. Propagator 

The propagator of antiprotons has a larger range than 
for positrons, and the effect of escape at the galactic 
boundaries must be considered. Modelling the galaxy 
as a cylinder of radius R ~ 20 kpc and of unknown 
height 2L, in which antiprotons diffuse with an energy- 
dependent coefficient K, the solution of the diffusion 
equation can be obtained as Bessel-Fourier expansion 
over the r and z variables. It turns out that the side 
boundaries (r = R) can be neglected in most situations, 
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FIG. 6: Distribution of the boost factors at 1 GeV (left panel) and 300 GeV (right panel) obtained with the Monte Carlo 
simulations and comparison to analytical expectations (see text for further details). The gaussian distribution discussed in the 
text is plotted in the left panel (long-dashed curve). 



as the z = ±L boundaries are closer to the solar neighbor- 
hood and influence more strongly the measured cosmic 
ray flux. As a consequence, the problem is essentially 
two-dimensional and the propagator from any source to 
the Earth can be expressed as a fonction or r and z only. 
Moreover, the sources (the mini-spikes) can be consid- 
ered as point-like as far as the propagation over galactic 
distances is considered, and special care has to be taken 
for the singular 1/r dependance of G(r,z). This point 
has been exposed in detail in [65], along with the cor- 
responding analytical expressions for G(r,z). The main 
difference with the positron case is that the energy losses 
due to radiative losses are much less important for an- 
tiprotons. If the diffusion is supposed to occur at con- 
stant energy, then the spatial distribution of antiprotons 
created in a given mini-spike is the same at every en- 
ergy, and the boost factor should be energy-independent 
(see below). There are actually at least two effects that 
depend on energy, namely spallation and galactic wind. 
Their influence will be discussed below, along with the 
results. We have considered the same median parameters 
as in the positron case. 



B. The Monte-Carlo simulation 

Following the same method as for positrons, wc have 
performed Monte-Carlo simulations of the IMBHs distri- 
bution, in order to estimate the expected antiproton flux, 
along with the associated variance. The procedure is ex- 
actly the same as before, and this section is devoted to 
the presentation of the results. 



C. Expected values and variance of the boost 
factor 



The average boost factor and its scatter are computed, 
both by direct evaluation of Eqs. (12) and (16) as well 
as by the Monte-Carlo simulation. Since energy losses 
are less important for antiprotons than for positrons, we 
find that when spallation and galactic wind are neglected, 
the boost factor does not depend on energy - see Fig. 7. 
When galactic wind is considered, propagation becomes 
energy-dependent at low energy, and so does the boost 
factor - see the dotted line in Fig. 7. The uncertainty 
band widens significantly below a few GeV. The impor- 
tance of the galactic wind can be estimated through the 
Peclet number defined as Pe = V C L/K, the effect of the 
wind being negligible for small values of Pc. The diffusion 
coefficient if is a growing function of energy, so propaga- 
tion is dominated by diffusion at high energy, whereas the 
wind has a sizeable effect at low energy. As the galactic 
wind is directed outward from the disk, it prevents an- 
tiprotons from reaching the Earth. As a result, the range 
of diffusion is lowered and the number of sources actually 
contributing to the signal is also lowered. This explains 
the larger variance at low energy. 

We see that the uncertainty on the antiproton flux is 
about one order of magnitude in the best case, i.e. with 
no energy-dependent processes. When both spallation 
and wind are considered, the prediction of the boost fac- 
tor becomes meaningless below ~ 20 GeV. Notice how- 
ever that when the thickness of the diffusive halo is larger, 
more mini-spikes contribute to the signal at the Earth 
and the variance decreases at all energies. 
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FIG. 7: Expected value and variance of the boost factor of 
the antiproton signal as a function of kinetic energy, in the 
case of a dark matter particle with mass m x — 1 TeV and for 
different p propagation configurations. 



V. A FEW REALISTIC MODELS FOR THE DM 
PARTICLE 

We have so far considered a generic case for the DM 
particle, with a mass of 1 TeV and a typical weak-scale 
annihilation cross section. This part is devoted to the es- 
timate of the positron and antiproton exotic fluxes within 
specific particle physics models. We choose four typical 
DM particle candidates in the frameworks of supersym- 
metry and extra-dimensions. 

All of these match the constraints from collider ex- 
periments and relic density. The smooth halo is mod- 
elled by a NFW profile with a Solar dark matter den- 
sity pq = 0.3 GeV/cm 3 . In the previous computations 
(positronic line, energy independent diffusion of antipro- 
tons), the energy distribution at the source g(Es) was 
not relevant since it cancelled in the ratio <j) r /<j) s . In this 
section, we make predictions for the flux and we per- 
form the convolution of the propagator with the injection 
spectrum g(Es). This is true in particular in the case of 
positrons, for which energy losses are more significant 
than for antiprotons. The injection spectrum depends 
on the final states of the DM particle annihilation pro- 
cess. In all the models we considered, these spectra have 
been computed with PYTHIA [67]. The annihilation cross 
section plays a significant role as regards the S factor of 
Eq. (7) as well as the annihilation volume £ through the 
dependence of the maximum density in the core, given 
by Eq. (19). 

In the framework of the Minimal Supersymmetric 
Standard Model (MSSM) with a gravity driven super- 
symmetry breaking (mSUGRA), we choose two specific 
sets of parameters leading to two typical DM candidates. 



The latter are a bino and a Higgsino (respectively referred 
to as B and H in the following), depending on the field 
content of the lightest supersymmetric particle. They 
differ mainly in their couplings to the Standard Model 
particles and thus have different annihilation final states. 
The public code SuSpect [fill] is used to derive the weak- 
scale parameters from the mSUGRA inputs by solving 
the renormalization group equations. The cross sections 
for the different final states have been determined with 
the micrOMEGAs [68, 69] package. 

The first Kaluza-Klein model considers warped extra- 
dimensions with a SO(10) GUT, as described in [70]. In 
that context the DM particle is a right handed Dirac 
neutrino. Its stability arises from the conservation of a 
Z 3 symmetry - hence its name LZP for Lightest Z Par- 
ticle. To fulfill the relic density constrains, we choose 
Mkk — 6 TeV and m LZ p = 50GcV respectively for 
the Kaluza-Klein scale and the LZP mass. The second 
Kaluza-Klein model assumes universal extra-dimensions 
with R- 1 = ITcV [ ]. Then, the DM particle is the 
lightest particle with odd K-parity, which happens to be 
the first Kaluza-Klein excitation of the B° hyperchargc 
gauge boson, referred to as the B^ (or LKP) with mass 
mLKP = 1 TeV. 

All the relevant particle physics parameters and the 
inferred values regarding the maximum density and the 
mean annihilation volume (£) arc summarized in Tab. I. 
Monte Carlo simulations are performed in the context of 
these four models. The LZP being a Dirac particle, one 
has 8 = 1/4 in that case, whereas 5=1/2 otherwise. 
The results are presented in Fig. 8 for positrons and in 
Fig. 9 for antiprotons. The absolute fluxes are featured 
together with the actual data points. The positron flux 
measurements are taken from HEAT-e* [72] and HEAT- 
pbar experiments [73] and two independent analyses of 
AMS01 data [74, 75]. For antiprotons, we borrowed the 
results from [77, 78, 79, 80]. 

We first discuss the positron case. In Fig. 8, the dashed 
line is the unboosted, and therefore guaranteed, exotic 
flux in the framework of the considered models. In each 
case, the yellow (grey) area represents the flux uncer- 
tainty with IMBH mini-spikes. For clarity, the bright 
yellow contours stand for the a/2 level whereas the 1 a 
contours are indicated by the pale yellow regions. In 
the four models we have considered, the positron flux 
at 10 GeV is in the range between ~ 6 x 10" 6 and 
~3x 10~ 5 cm" 2 s" 1 GeV" 1 sr" 1 . The exotic fluxes gen- 
erated by IMBH mini-spikes arc fairly insensitive to the 
annihilation cross section, as a consequence of Eq. (22). 
We unexpectedly find that the flux is not very sensitive to 
the mass. This is because the injection spectrum g(Es) 
softens the m x dependence of the annihilation rate in- 
ferred from that relation. In addition, the variance of the 
flux is significantly smaller in the LKP case. Indeed, the 
injection spectrum g(E$) extends up to 1 TeV and has an 
important high energy contribution. As a consequence, 
the sensitivity sphere below ~10 GeV has a wider spread 
than for the three other DM candidates. The LKP flux 
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FIG. 8: Exotic positron fluxes in four particle physics models (top : supersymmetric DM, bottom : Kaluza-Klein DM) 
within IMBHs mini-spike scenario, and comparison to actual measurements. The bright and the pale yellow areas respectively 
correspond to contours at the a/2 and 1 a level. 



is less sensitive to the shotnoisc associated with the ran- 
dom realization of the IMBH population. The other cases 
feature a higher dispersion, mainly due to the fact that 
the DM particle mass is lower and the sensitivity sphere 
smaller. 

We now discuss the antiproton case. The results are 
displayed in Fig. 9, along with the standard flux, for the 
same set of propagation parameters as for the positrons 
(see Sec. III). In the case of a smooth halo, the exotic 
contribution is always lower than the flux due to stan- 
dard processes. When IMBH mini-spikes are considered, 
the antiproton flux at 10 GeV is in the range between 
-3 xl0~ 2 and -4 x lO" 1 irT 2 s" 1 GeV -1 sr" 1 in the four 
models we have considered. We observe the same deli- 
cate interplay between the mass dependence of the an- 
nihilation rate and the energy behaviour of the injection 
spectrum g(E$) - though the later is less dependent on 



the final state species than for positrons. As before, the 
antiproton flux is weakly dependent on the annihilation 
cross section. This is particularly obvious in the two up- 
per panels, between which ((Tannf) varies by a factor ~7. 
Below ^20 GeV, the antiproton primary fluxes have a 
very large scatter, as explained in Sec. IV C. 



VI. DISCUSSION AND PERSPECTIVES 

We have considered a scenario in which the formation 
of IMBHs generates dramatic enhancements of the DM 
density in their vicinity. We have computed the expected 
number of unmerged IMBHs in the Galactic halo, their 
spatial distribution and the properties of their associated 
mini-spikes. We found high values for the positron and 
antiproton boost factors, e.g. of order 10 4 for a fiducial 



11 



: = 140 GeV 

■ smooth contribution 

- bg from Maurin et al 

- smooth+mini-spikes+bg 
modulated at <|> = 500 MV 



BESS 95+97 
BESS 98 
BESS 99 
BESS 00 



WIZARD-MASS 91 
CAPRICE 94 
CAPRICE 98 
AMS 98 




T p [GeV] 



> 



e 



io- 6 r 
iff 7 — 
iff 8 =- 



biggsino: m = 108 GeV 

smooth contribution 

bg from Maurin et al 

smooth+mini-spikes+bg 

modulated at <|> = 500 MV 



10 1 



BESS 95+97 
BESS 98 
BESS 99 
BESS 00 



WIZARD-MASS 91 
CAPRICE 94 
CAPRICE 98 
AMS 98 




T p [GeV] 



LZP:m ( = 50 GeV 

smooth contribution 

bg from Maurin et al 

smooth+mini-spikes+bg 

modulated at <|> = 500 MV 



BESS 95+97 
BESS 98 
BESS 99 
BESS 00 



WIZARD-MASS 91 
CAPRICE 94 
CAPRICE 98 
AMS 98 




> s 

o 

IN 

E 1 



O 10 " : 
■o 



!.,= 1000 GeV 
■ smooth conlribution 

- bg from Maurin et al 

- smooth+mini-spikes+bg 
modulated at <|> = 500 MV 



WIZARD-MASS 91 
CAPRICE 94 
CAPRICE 98 
AMS 98 




T B [GeV] 



T B [GeV] 



FIG. 9: Antiproton fluxes as a function of kinetic energy in the case of a bino dark matter particle with mass m x = 140 GeV. 
The fluxes were computed without (upper panel) or with (lower panel) wind (with a value V c = 12 km/s), without (left) or 
with (right) spallation. 



DM particle with a 1 TeV mass. For an average mini- 
spike population, the flux of primary positrons and an- 
tiprotons turns out to be one to two orders of magnitude 
larger than the measured flux, and is fairly insensitive 
to the specific properties of the DM candidate. In par- 
ticular, the annihilation cross section has little influence 
because of the presence of an annihilation plateau in the 
cores of mini-spikes. This interesting feature allows con- 
figurations in the particle physics parameter space which 
are usually disregarded - because of their very low cross 
sections - to become testable. As the average flux signal 
exceeds the data points, we could naively conclude that 
the mini-spike scenario is already ruled out by observa- 
tion. However, because of the small number of objects, 



the variance associated to the positron and antiproton 
signal is enormous and no definitive conclusion can be 
reached. We leave for further investigation the quanti- 
tative estimate of the confidence level at which the the 
mini-spike scenario agrees with measurements. 

Notice finally that this variance is small at low energy 
for positrons, and on the contrary at high energy for an- 
tiprotons. Should a cut-off above some specific energy 
E c be detected in the positron flux, the natural interpre- 
tation would lead to a DM particle mass m x ~ E c . The 
possible existence of substructures leaves room to more 
subtle explanations. In particular, if an excess in the an- 
tiproton flux were found above E C1 we would conclude 
that we live in a particular Milky Way realization where 
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substructures are far from the Earth. The farther those 
substructures, the smaller E c with respect to m x , hence 
the possibility of a misinterpretation of the data. This il- 
lustrates the fact that combining various indirect signals 
could give information about the spatial distribution of 
DM. 
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B 


H 


LZP 


LKP 


m 


140 GeV 


108 GeV 


50 GeV 


1 TeV 


(7v (xlO -26 em's" 1 ) 


0.26 


1.9 


2.04 


1.7 


Final States 


66 91% 

r+T- 9% 


W+W- 90% 
Z°Z° 10% 


qq 74% 
vv 17% 
(.+£- 2.88% (x3) 


qq up 11% (x3) 
qcjdown 1% (x3) 
vv 4% 

e+r 20% (x3) 


p max (xlO 17 M kpc' 3 ) 


47.1 


4.97 


2.15 


51.5 


(0 (xlO 5 kpc 3 ) 


46.6 


9.35 


5.14 


49.6 



TABLE I: Relevant parameters for the particle dark matter models discussed in Sec. V. 



